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SCHRODINGER OPERATORS WITH SINGULAR PERIODIC 

POTENTIALS 
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Abstract. By using quasi-derivatives, we develop a Fourier method for 
studying the spectral properties of one dimensional Schrodinger operators 
with periodic singular potentials. 



1. Introduction 

Our goal in this paper is to develop a Fourier method for studying the spectral 
properties (in particular, spectral gap asymptotics) of the Schrodinger operator 

(1.1) L{v)y = —y" + v{x)y, x G M, 
where v is a singular potential such that 

(1.2) v{x) = v{x + tt), ve Hi-^l{R). 

In the case where the potential is a real L^([0, 7r])-function, it is well known 
by the Floquet-Lyapunov theory (see [5l HHJ EOj |32] ) , that the spectrum of L is 
absolutely continuous and has a band-gap structure, i.e., it is a union of closed 
intervals separated by spectral gaps 

(-oo,Ao), (A5",A+), (A2,A+),--- ,(A^,A+),--- . 

The points (A^) are defined by the spectra of (jl.ip considered on the inter- 
val [0,7r], respectively, with periodic (for even n) and anti-periodic (for odd n) 
boundary conditions (be) : 

(a) periodic Per+ : yi^n) = y{0), y'(7r) = y'{0); 

(b) antiperiodic Per~ : y{Tr) = — y(0), y'(vr) = — y'(0); 

So, one may consider the appropriate bases in L^([0, vr]), which leads to a 
transformation of the periodic or anti-periodic Hill-Schrodinger operator into an 
operator acting in an ^^-sequence space. This makes possible to develop a Fourier 
method for investigation of spectra, and especially, spectral gap asymptotics 
(see [m [H], where the method has been used to estimate the gap asymptotics 
in terms of potential smoothness). Our papers [21 [3] (see also the survey [4]) 
give further development of that approach and provide a detailed analysis of 
(and extensive bibliography on) the intimate relationship between the smoothness 
of the potential v and the decay rate of the corresponding spectral gaps (and 
deviations of Dirichlet eigenvalues) under the assumption v G L^([0,7r]). 

But now singular potentials v G bring a lot of new technical problems 
even in the framework of the same basic scheme as in [4] . 
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First of them is to give proper understanding of the boundary conditions (a) 
and (b) or their broader interpretation and careful definition of the corresponding 
operators and their domains. This is done by using quasi-derivatives. To a great 
extend we follow the approach suggested (in the context of second order o.d.e.) 
and developed by A. Savchuk and A. Shkalikov [25l[27| (see also [26| [28l [29] ) and 
R. Hryniv and Ya. Mykytyuk [8] (see also [9]-[l3]). For specific potentials see 
W. N. Everitt and A. Zettl \6',^. 

E. Korotyaev [171 [TSj follows a different approach but it works only in the case 
of a real potential v. 

It is known (e.g., see Remark 2.3, or Proposition [1] below) that every 
TT-periodic potential v G H^^ (M) has the form 

V = C + Q' , where C = const, Q is vr — periodic, Q G L^^^(M). 

Therefore, one may introduce the "quasi-derivative" u = y' — Qy and replace 
the distribution equation —y" + vy = by the following system of two linear 
equations with coefficients in Lj^^{M.) 

(1.3) y' = Qy + u, v! = [C - Q^)y - Qu. 

By the Existence-Uniqueness theorem for systems of linear o.d.e. with L^^^^(M)- 
coefhcients (e.g., see [Il[22]), the Cauchy initial value problem for the system ()1.3p 
has, for each pair of numbers (a, 6), a unique solution (y,u) such that y(0) = 
a, tt(0) = h. 

Moreover, following A. Savchuk and A. Shkalikov }25|, I27j. one may consider 
various boundary value problems on the interval [0, vr]). In particular, let us 
consider the periodic or anti-periodic boundary conditions Per^, where 

(a*) Per+ : y(7r) = y(0), {y' - Qy) (vr) = {y' - Qy) (0). 

{b*) Per- : ^(vr) = -y(0), (y' - Qy) (vr) = - {y' - Qy) (0). 

R. Hryniv and Ya. Mykytyuk [8j used also the system (jl.3p in order to give 
complete analysis of the spectra of the Schrodinger operator with real-valued 
periodic ^-potentials. They showed, that as in the case of periodic ^^^^(M)- 
potentials, the Floquet theory for the system (II. 3p could be used to explain that if 
V is real-valued, then L(v) is a self-adjoint operator having absolutely continuous 
spectrum with band-gap structure, and the spectral gaps are determined by the 
spectra of the corresponding Hill-Schrodinger operators Lp^^± defined in the 
appropriate domains of L^([0, 7r])-functions, and considered, respectively, with 
the boundary conditions (a*) and (6*). 

In Section 2 we use the same quasi-derivative approach to define the domains 
of the operators L{v) for complex-valued potentials v, and to explain how their 
spectra are described in terms of the corresponding Lyapunov function. From 
a technical point of view, our approach is different from the approach of R. 
Hryniv and Ya. Mykytyuk [8j: they consider only the self-adjoint case and use 
a quadratic form to define the domain of L{v), while we consider the non-self- 
adjoint well and use the Floquet theory and the resolvent method (see 
Lemma [3] and Theorem [4]). 



FOURIER METHOD 



3 



Sections 3 and 4 contains the core results of this paper. In Section 3 we 
define and study the operators Lper± which arise when considering the Hill- 
Schrodinger operator L{v) with the adjusted boundary conditions (a*) and (b*). 
We meticulously explain what is the Fourier representation of these operatorq^ 
in Proposition [10] and Theorem II 1[ 

In Section 4 we use the same approach as in Section 3 to define and study 
the Hill-Schrodinger operator L£)ir{v) with Dirichlet boundary conditions Dir : 
y(0) = y(vr) = 0. Our main result there is Theorem [16] which gives the Fourier 
representation of the operator L£)ir{v). 

In Section 5 we use the Fourier representations of the operators Lp^^± and 
L£)ir to study the localization of their spectra (see Theorem 15. ip . Of course, 
Theorem 15.11 gives also a rough asymptotics of the eigenvalues , A~ , of 
these operators. But we are interested to find the asymptotics of spectral gaps 
7n = — A~ in the self-adjoint case, or the asymptotics of both 7„ and the 
deviations — (A^ + A~)/2 in the non-self-adjoint case, etc. Our results in that 
direction are presented without proofs in Section 6. 

Acknowledgment. The authors thank Professors Rostyslav Hryniv, Andrei 
Shkalikov and Vadim Tkachenko for very useful discussions of many questions of 
spectral analysis of differential operators, both related and unrelated to the main 
topics of this paper. 

2. Preliminary results 

1. The operator 11.11 has a second term vy with v E (jl.2p . First of all, let 
us specify the structure of periodic functions and distributions in Hl^^{M.) and 

The Sobolev space H^^J^) is defined as the space of functions f{x) G L^^^(M) 
which are absolutely continuous and have their derivatives f'{x) G L^^^(]R). 
Therefore, for every T > 0, 

(2.1) ||/||1t = 1^ (1/(^)1' + |/'(x)H <oo. 

Let D(M) be the space of all C°°-functions on M with compact support, and 
let P([-r,r]) be the subset of all ip E P(M) with suppip C [-T,T]. 
By definition, H~l{m.) is the space of distributions on M such that 

(2.2) Vr>0 3C7(r): \{v,ip)\<C{T)M\i,T y^ev{[-T,T]). 

^Maybe it is worth to mention that T. Kappeler and C. Mohr [16] analyze "periodic and 
Dirichlet eigenvalues of Schrodinger operators with singular potential" but they never tell how 
these operators (or boundary conditions) are defined on the interval, i.e., in a Hilbert space 
L2([0,7r]). At some point they jump without any justification or explanation into weighted 
£^-sequence spaces (an analog of Sobolev spaces ff") and consider the same sequence space 
operators we are used to in the regular case, i.e., if u £ I/pg^(K). But without formulating which 
Sturm-Liouville problem is considered, what are the corresponding boundary conditions, what 
is the domain of the operator, etc., it is not possible to pass from a non-defined differential 
operator to its Fourier representation. 
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Of course, since 

T r-T 

\ip{x)\^dx <T^ \ip'{x)\'^dx, 

-T J~T 

the condition (12. 2p is equivalent to 

(2.3) Vr>0 3C7(r): \{v,V^)\<CiT)\\ip'\\L2^[_T,T]) V99 G p([-r, T]). 
Set 

(2.4) Pi(M) = V: ipeV{R)}, Vi{[-T,T]) = y : ^eV{[-T,T])} 
and consider the hnear functional 

(2.5) q{ip'):=-{v,ip), ^' £ Vi( 



In view of (|2.3p . for each T > 0, g(-) is a continuous linear functional defined in 
the space Di{[—T,T]) C L^([— T, T]). By Riesz Representation Theorem there 
exists a function Qt{x) G L^([— T, T]) such that 

(2.6) q{'p') = j QT{xW{x)dx V(^eP([-r,r]). 

The function Qt is uniquely determined up to an additive constant because in 
L^([— r, T]) only constants are orthogonal to 'Di{[—T,T]). Therefore, one can 
readily see that there is a function Q{x) G Lij^JW) such that 

/oo 
Q{x)ip'{x)dx V9? G P(1R), 
-00 

where the function Q is uniquely determined up to an additive constant. Thus, 
we have 

{v,ip) = -q{(p') = -{Q,(p') = {Q',ip), 

i.e., 

(2.7) V = Q'. 

A distribution v G H^^ (M) is called periodic of period vr if 

(2.8) {v,ip{x)) = {v,ip{x-TT)) V(^gP(IR). 

L. Schwartz ^30j gave an equivalent definition of a periodic of period vr distri- 
bution in the following way: Let 

u! : W ^ = M/vrZ, uj{x) = x mod vr. 

A distribution F G P'(M) is periodic if, for some / G (C°°(S'"^))' , we have 

F{x) = f{u;{x)), i.e., (<^, F) = (cD, /), 

where 

$ = (/?(x — /cvr). 
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Now, if V is periodic and Q G L^^^(M) is chosen so that (|2.7p holds, we have 
IHD 

pcio poo fOO 

Q{x + Tr)ip' {x)dx = / Q{x)(p' {x — tt) = / Q{x)(p' {x)dx, 



I.e., 

/oo 
[Q{x + TT)-Q{x)]ip'{x)dx = V(^eP 
-oo 

Thus, there exists a constant c such that 

Q(x + vr) — (5(x) = c a.e. 

Consider the function 

c 

= x; 

IT 

then we have Q{x + tt) = Q{x) a.e., so Q is vr-periodic, and 

v = Q' + -. 

vr 

Let 

(2.9) Q(x) = (?(m)e™- 

me2Z 

be the Fourier series expansion of the function Q G L^([0,7r]). Set 

(2.10) y(0) = -, V{m) = imq{m) for m / 0. 

vr 

Ah this leads to the following statement. 

Proposition 1. Every t: -periodic distribution v G ff^~^(M) has the form 

(2.11) v = C + Q', QgLL(M), Q(x + ^)=-Q(x) 

1 r 

(2.12) ^(0) = - / Q{x)dx = 0, 

71" Jo 

and can he written as a converging in i?^^](]R) Fourier series 

(2.13) v=Y^ V{m)e''^^ 

with 

(2.14) V{0) = C, V{m) = imq{m) for m / 0, 
where q{m) are the Fourier coefficients ofQ. Of course, 

(2.15) IIQIll,,,,.,, = E 
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Remark. R. Hryniv and Ya. Mykytyuk [8j, (see Theorem 3.1 and Remark 2.3) 
give a more general claim about the structure of uniformly bounded if^~^(M)- 
distributions. 

2. In view of (|2.2p . each distribution v G H^^{M.) could be considered as a 
linear functional on the space Hlg{M.) of functions in H^^^(K) with compact sup- 
port. Therefore, if u G H^^CR) and y G Hj-^^CR), then the differential expression 
£{y) = —y" + V ■ y is well-defined by 

{-y" + vy,ip) = {y', ip') + {v,y- Lp) 

as a distribution in //^^^(M). This observation suggests to consider the Schrodinger 
operator —d^/dx^ + v in the domain 

(2.16) D{L{v)) = {y^Hl,{R)nL\R) : -y" + v ■ y ^ L\R)] . 

Moreover, suppose v = C + Q' , where C is a constant and Q is a vr-periodic 
function such that 

(2.17) Q£L\[0,tt]), q{0) = - [ Q{x)dx = 0. 

Jo 

Then the differential expression £{y) = —y" + vy can be written in the form 

(2.18) £(y)=_(y'_Qy)'_Qy' + C7y. 

Notice that 

iiy) = -{y- Qy)' -Qy' + Cy = fG 

if and only if 

u = y'-QyG Wli,M) 
and the pair {y, u) satisfies the system of differential equations 



(2.19) 



y' = Qy + u, 

u' = {C- Q^)y -Qu + f. 



Consider the corresponding homogeneous system 

y' = Qy + u, 
u' = {C- Q'^)y - Qu. 



(2.20) 
with initial data 

(2.21) y(0) = a, u{0) = b. 

Since the coefficients 1,Q,C — of the system (j2.20p are in Lj^^{R), the stan- 
dard existence-uniqueness theorem for linear systems of equations with Lj^^(K)- 
coefficients (e.g., see M. Naimark [22], Sect. 16, or F. Atkinson [1]) guarantees 
that for any pair of numbers (a, b) the system (|2.20p has a unique solution {y, u) 
with y,u£ Wli^^{R) such that d^Tli holds. 

On the other hand, the coefficients of the system (12.201) are vr-periodic, so one 
may apply the classical Floquet theory. 
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Let and {1/2, U2) be the solutions of ()2.20p which satisfy yi(0) = 

l,Mi(0) = and y2(0) = 0,U2(0) = 1. By the Caley-Hamilton theorem the 
Wronskian 



det 



f yi{x) y2{x)\ _ ^ 

\Ui{x) U2{x)J 



because the trace of the coefficient matrix of the system (|2.20p is zero. 

If {y{x), u{x)) is a solution of ()2.20p with initial data (a, b), then {y{x+TT),u{x+ 
tt)) is a solution also, correspondingly with initial data 

u(jt)J \hj \ui[-k) U2{7r) 

Consider the characteristic equation of the monodromy matrix M : 

(2.22) p'^-Ap+l = 0, A = yi(7r) +U2(7r). 

Each root p of the characteristic equation (12.22P gives a rise of a special solution 
{ip{x),i;{x)) of (f2:20]l such that 

(2.23) ip{x + tt) = p ■ ip{x), '4){x + vr) = p • ^(x). 

Since the product of the roots of (j2.22p equals 1, the roots have the form 

(2.24) p^ = e^^^, T = a + if3, 

where /? S [0, 2] and a = if the roots are on the unit circle or a > otherwise. 

In the case where the equation (I2.22P has two distinct roots, let {ip^,ip^) be 
special solutions of (12.201) that correspond to the roots (|2.24p . i.e., 

(93=^(x + 7r),V'=^(x + 7r)) = p^ ■ (cp^ (x) , tp^ (x)) . 

Then one can readily see that the functions 

(x) = e^"^v^± (x) , ^± (x) = e^™' (x) 

are vr-periodic, and we have 

(2.25) (^±(x) = e±™(^±(x), ?/;±(x) = e±™V^±(x). 

Consider the case where ()2.22p has a double root p = ±1. If its geometric mul- 
tiplicity equals 2 (i.e., the matrix M has two linearly independent eigenvectors), 
then the equation (j2.20p has, respectively, two linearly independent solutions 
(99^,^^) which are periodic if p = 1 or anti-periodic ii p = —1. 

Otherwise, (if M is a Jordan matrix), there are two linearly independent vec- 
tors (^+^ ) such that 
(2.26) 

" ) - " ■ {:-) - " (:-) + ) ■ ^ - ^ »■ 

Let {(p^,i/j'^) be the corresponding solutions of (j2.20p . Then we have 
(2.27) 

(P^{x + 'it)\ j ip^(x)\ jip^(x + 7r)\ j ip^{x)\ (ip'^{x) 



iP+ix + T:)] ~'"\iP+{x)J ' \ip-{x + TT)J ~^U"(x); ^P^\i;+{x) 
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Now, one can easily see that the functions (p and "tp given by 

•f'~{x)\ ^ / v~{x)\ <^^{xy 

are 7r-periodic (if p = 1) or anti-periodic (if p = —1). Therefore, the solution 

, _ ; ( can be written in the form 
{x)J 

^ ' {i^-ix)J ~ \^P~ix)J + TT [i^+ix)^ 

i.e., it is a linear combination of periodic (if /> = 1), or anti-periodic (if p = —1) 
functions with coefficients 1 and kx/tt. 

The following lemma shows how the properties of the solutions of ()2.19p and 
()2.20p depend on the roots of the characteristic equation (j2.22p . 

Lemma 2. (a) The homogeneous system i2. 20\) has no nonzero solution {y,u) 
with y G L^(M). Moreover, if the roots of the characteristic equation i2. 22\) lie 
on the unit circle, i.e., a = ^ in the representation ( 2.24% then \2.20\) has no 
nonzero solution {y,u) with y G L^((— oo,0]) or y ^ L^([0, +oo)). 

(b) If a = in the representation 1^2.24% then there are functions f S L^(M) 
such that the corresponding non-homogeneous system ^2.19\) has no solution (y, u) 
with y £ L2(M). 

(c) If the roots of the characteristic equation 12. 2^) lie outside the unit circle, 
i.e., a > in the representation \2. 24\ ), then the non-homogeneous system i2.19\) 
has, for each f £ L^(IR), a unique solution {y,u) = (i?i(/), i?2(/)) such that 
Ri is a linear continuous operator from L^(]R) into W2^(M), and R2 is a linear 
continuous operator in L^(M) with a range in W^i^^iM.). 

Proof, (a) In view of the above discussion (see the text from (12.221) to (12.28P ). 
if the characteristic equation (j2.22p has two distinct roots p = e^'^'^, then each 
solution (y, u) of the homogeneous system (j2.20p is a linear combination of the 
special solutions, so 

y{x) = C+e^^(^+(x) + C-e-^^'(^-(x), 

where (p'^ and ip~ are 7r-periodic functions in H^. 

In the case where the real part of r is strictly positive, i.e., r = a + ifi with a > 
0, one can readily see that e™(^+(a;) I/^([0,oo)) but e^^(^"'"(x) G L^((— oo,0]), 
while e-™(^-(a;) G L2qo,oo)) but e-™(^-(x) L2((-oo, 0])). Therefore, if y ^ 
we have y L^(M). 

Next we consider the case where t = i(3 with /? 7^ 0, 1. The Fourier series of 
the functions <f~^{x) and (p~{x) 

fce2Z fce2Z 
converge uniformly in M because (p~^,ip^ G H^. Therefore, we have 

k&2Z k&2Z 
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where the series on the right converge uniformly on M. If /? is a rational number, 
then y is a periodic function, so y ^ cxd, 0]) and y -^^^([0, oo)). 

If (3 is an irrational number, then 

1 

(2.29) lim - / y(2;)e-*('=^^)^dx = C^^^ VA: G 2Z. 

T-*oo T Jq 

On the other hand, if y G -^^([0) oo)), then the Cauchy inequality implies 

1 / /.r \ 1/2 / \ 1/2 II II 



But, in view of (j2.29p . this is impossible if y 7^ 0. Thus y L^([0, cxd)). In a 
similar way, one can see that y L2((_oo^o]). 

Finally, if the characteristic equation (j2.22p has a double root p = ±1, then 
either every solution (y, u) of (j2.20p is periodic or anti-periodic, and so y 
L2([0,oo) and y ^ L2((-oo,0]), or it is a linear combination of some special 
solutions (see (|2.28p . and the preceding discussion), so we have 

y(x) = C+ip+ix) + C-^- + C-—^+{x), 

vr 

where the functions ip~^ and (f~ are periodic or anti-periodic. Now one can easily 
see that y L^([0, 00) and y L^((— oo,0]), which completes the proof of (a). 

(b) Let {ip^,xp^) be special solutions of (|2.20p that correspond to the roots 
(j2.24p as above. We may assume without loss of generalities that the Wron- 
skian of the solutions {ip^,^p^) and {ip~ ,xp~) equals 1 because these solutions are 
determined up to constant multipliers. 

The standard method of variation of constants leads to the following solution 
(y, n) of the non-homogeneous system (j2.19p : 

(2.30) y = {x)ip'^ (x) + {x)ip^ (x), u = {x)Tp^ (x) + v~ {x)7p~ (x), 
where and v~ satisfy 

dv^ , dv^ dv'^ , dv~ ^ 

2.31 — •v'+ + :r •¥'"=0' T- -"f" +T- •V'" = /, 

dx dx dx dx 

so 

(2.32) v^{x) = - ip"{t)f{t)dt + C^, v-{x)= ip+{t)f{t)dt + C-. 
Jo Jo 

Assume that the characteristic equation (j2.22p has roots of the form p = 

gi/37r^ P g 2). Take any function / G L^(IR) with compact support, say supp/ C 

(0, T). By (j2.30p and (|2.32p . if (y, u) is a solution of the non-homogeneous system 

()2.19p . then the restriction of {y,u) on the intervals (— oo,0) and [T, cx)) is a 

solution of the homogeneous system l^20h . So, by (a), if y G L^{R) then y = 

on the intervals (— oo,0) and [T, 00). This may happen if only if the constants 

in (|2.32p are zeros, and we have 

r it)fit)dt = 0, r ip+{t)f{t)dt = 0. 

Jo Jo 
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Hence, if / is not orthogonal to the functions 99^, on the interval [0, T], then 
the non-homogeneous system ()2.19p has no solution {y,u) with y € L'^(M). This 
completes the proof of (b). 

(c) Now we consider the case where the characteristic equation (j2.22p has 
roots of the form (|2.24p with a > 0. Let {ip^,tl^^) be the corresponding special 
solutions. By ([230]) . for each / G ^^(M), the non-homogeneous system ()2.19p 
has a solution of the form (y,n) = i?2(/); where 

(2.33) 

Riif) = {x)lp^ {x) + {x)(p^ {x), R2{f) = {x)'ip~^ (x) + v~ {x)'ip~ (x) , 

and (j2.3ip holds. In order to have a solution that vanishes at ±00 we set (taking 
into account (j2.25p l 

(2.34) v+{x)= e~^^^~{t)f{t)dt, v' {x) = e^^^+ {t)f{t)dt. 

J X J — 00 

Let C± = max{|(^±(x)| : x G [0,7r]}. By (^M), we have 



(2.35) \if^{x)\ <C±-e 

Therefore, by the Cauchy inequality, we get 

2 



\v+(x)\' < Ci 



00 \ / fOO 



<Ci{ e-'^'dt • / e-'^'\f{t)\'dt 



so 



fi2 

(2.36) l^'^(^)l^ < — e-"^' / e-"'|/(OI'dt. 

a 

Thus, by (l235l) . 



00 poo poo 

\v+{x)r\ip+{x)rdx<—^ e-'''\fit)\'^dtdx 



^'^ \ ip' (x)\ ax s — 
00 ^ J —00 



/^<2 r^2 foo / rt \ r^2 

<^ \m\H e"(-*)dx ) = ^ 
In an analogous way one may prove that 

\v- {x)f (x)f dx <^^yf^,. 



a 

2 



In view of (j2.30p , these estimates prove that R\ is a continuous operator in L 
Next we estimate the L^(M)-norm of y' = -^Riif)- In view of (I2.3ip . we have 

// \ 4-/ \ dio^ , . , dip^ , . 

y (x) = V (x) ■ — ix) + V [X) ■ — — (x). 

dx dx 



By (ESSD, 



v+{x) ■ ^^(x) = av+{x)ip+ + 7;+(x)e"^^ 
dx dx 
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Since the L^(M)-norm of v^{x)ip^ has been esthnated above, we need to estimate 
only the L^(M)~norm of {x)e'^^ d(f~^ / dx . By (I2.36p . we have 

/oo /~^2 POO POO 

\v+{x)e'^''d^+/dx\^dx<— e"^ / e 

-oo ^ J ~oo Jx 



\f{t)\^dtdx 



a 



\f{t)? 



|#+/(ix|^e°(^~*)dx ) dt. 



Firstly, we estimate the integral in the parentheses. Notice that the function 
d^p^ldx (and therefore, d(p^/dx ) are in the space L2([0,7r]) due to the first 
equation in (|2.2U|) . Therefore, 



(2.37) 
We have 



dip 
dx 



(x) 



dx < oo. 



Q J-(n+l)7r 



dip' 
dx 



e'^^di 



< 



n=0 



-<(l + a7r)^ 
1 — exp(— avrj avr 



Thus, 



v+{x)e'^''dip+ /dx\ dx < (1 + avr) 



In an analogous way it follows that 

.-(x)e"-^"' 
dx 



dx < {1 + avr 



so the operator Ri act continuously from L^(]R) into the space 

The proof of the fact that the operator R2 is continuous in L^(M) is omitted 
because essentially it is the same (we only replace 99^ with ijj^ in the proof that 
i?i is a continuous operator in L^(]R)). □ 

We need also the following lemma. 

Lemma 3. Let H he a Hilhert space with product (•, •), and let 

A : D{A) ^ H, B : D{B) H 

be (unbounded) linear operators with domains D(A) and D{B), such that 

(2.38) {Af, g) = (/, Bg) for f G D{A), g G D{B). 

If there is a X £ C such that the operators A — X and B — X are surjective, then 
(i) D[A) and D{B) are dense in H; 

(a) A* = B and B* = A, where A* and B* are, respectively, the adjoint 
operators of A and B. 
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Proof. We need to explain only that D{A) is dense in H and A* = B because 
one can replace the roles of A and B. 

To prove that D{A) is dense in H, we need to show that if h is orthogonal to 
D{A) then h = 0. Let 

{f,h) = o yfeD{A). 

Since the operator B — X is surjective, there is g & D{B) such that h = {B — X)g. 
Therefore, by (|2.38p . we have 

= (/, h) = if, (B - X)g) = {{A - A)/, g) V/ G D{A), 

which yields g = because the range of A — A is i^. Thus, h = {B — X)g = 0. 
Hence (i) holds. 

Next we prove (ii). If g* G Dom{A*), then we have 

(2.39) {A-X)f,g*) = if,w) yfGD{A), 

where w = {A* — X)g* . Since the operator B — X is surjective, there is 5 G D[B) 
such that w = {B — X)g. Therefore, by (I2.38|) and ()2.39p . we have 

{{A - X)f,g*) = if, iB-X)g) = {{A - X)f,g) V/ G D{A), 

which implies that g* = g (because the range of A — A is equal to H) and 
{A* - X)g* = {B - X)g*, i.e., A*g* = Bg*. This completes the proof of (ii). 

□ 

Consider the Schrodinger operator with a spectral parameter 

L{v) - A = -(f/dx^ + {v-X), AG C. 

In view of the formula (j2.11l in Proposition [H we may assume without loss of 
generality that 

(2.40) C7 = 0, v = Q', 

because a change of C results in a shift of the spectral parameter A. 
Replacing C by —A in the homogeneous system (j2.20p . we get 



(2.41) 



y' = Qy + u, 

u' = {-X-Q^)y-Qu. 



Let {yi{x; X),ui{x; X)) and (2/2(2;; A), M2 (2;; A)) be the solutions of (j2.4ip which 
satisfy the initial conditions yi(0; A) = 1, ni(0; A) = and 2/2(0; A) = 0, U2{0; A) = 
1. Since these solutions depend analytically on A G C, the Lyapunov function, or 
Hill discriminant, 

(2.42) A(Q,A) =yi(7r;A) + U2(vr;A) 

is an entire function. Taking the conjugates of the equation in (I2.41|] . one can 
easily see that 



(2.43) A(Q,A) = A(Q,A). 
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Remark. A. Savchuk and A. Shkalikov gave asymptotic analysis of the func- 
tions yj(7r,A) and Uj{7r,X), j = 1,2. In particular, it follows from Formula (1.5) 
of Lemma 1.4 in [27] that, with = A, 
(2.44) 

yi(7r,A) = cos(7rz)-|-o(l), y2(7r,A) = -[sin(7rz)-|-o(l)], n2(7r, A) = cos 7rz-|-o(l), 
and therefore, 

(2.45) A(Q,A) = 2cos7rz + o(l), z"^ = X, 
inside any parabola 

(2.46) Pa = {XeC: \Im z\ < a}. 

In the regular case v £ L^([0,7r]) these asymptotics of the fundamental solutions 
and the Lyapunov function A of the Hill-Schrodinger operator could be found 
in p. 32, Formula (1.3.11), or pp. 252-253, Formulae (3.4.23'), (3.4.26). 

Consider the operator L{v), in the domain 
(2.47) 

D{L{v)) = {y e H\R) : y' - Qy e (R) H Wli,,{R) , £Q{y) e L^iR)} , 
defined by 

(2.48) L{v)y = £Q{y), with igiy) = -{y' - Qy)' - Qy' , 
where v and Q are as in Proposition [TJ 

Theorem 4. Let v G Hj^^{R)) be ir-periodic. Then 

(a) the domain D{L{v)) is dense in L^(M); 

(b) the operator L{v) is closed, and its conjugate operator is 

(2.49) {Liv)r = Hv); 

(In particular, if v is real-valued, then the operator L{v) is self -adjoint.) 

(c) the spectrum Sp{L[v)) of the operator L{v) is continuous, and moreover, 

(2.50) Sp{L{v)) = {A G C I 30 G [0, 2tt) : A(A) = 2 cos 6]. 

Remark. In the case of L^-potential v this result is known (see Rofe-Beketov 
[231 [Ml and V. Tkachenko [3l]). 

Proof. Firstly, we show that the operators L{v) and L{v) are formally adjoint, 
i.e., 

(2.51) {L{v)y, h) = (/, L{v)h) if y G D{L{v)), h G D{L{v)). 

Since y' — Qy and h are continuous L^(]R)-functions, their product is a continuous 
-L-'^(M)-function, so we have 

liminf \ {y' — Qy)h\ (x) = 0. 

Therefore, there exist two sequences of real numbers c„ ^ — oo and (i„ — > oo 
such that 

{{y - QyW) (cn) 0, ((y' - Qy)h) (dn) ^0 as n ^ oo. 
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Now, we have 

{L{v)y,h)= I i.Q{y)hdx= lim / {-{y - Qy)'h - Qy'h) dx 

dn rdn rdn 



lim 

n— >oo 



"m ran ran ' 

{y' - Qy)h + {y' - Qy)h'dx - / Qy'hdx 

J Cji Crj, I 



/OO 
{y'W - Qyhl - Qy'h) dx. 
-oo 

The same argument shows that 

[y'h' - Qyh' - Qy'h) dx = {y, L{v)h) , 

which completes the proof of (j2.5ip . 

If the roots of the characteristic equation — A((5, X)p+ 1 = lie on the unit 
circle {e*^, 6 S [0, 27r)}, then they are of the form e^*^, so we have 

(2.52) A(Q,A) = e^^ + e"^^ = 2cos0. 

Therefore, if A((5, A) [—2,2], then the roots of the characteristic equation lie 
outside of the unit circle {e^^ ,6 G [0,27r)}. If so, by part (c) of Lemma [21 the 
operator L{v) — A maps bijectively D{L{v)) onto L^(M), and its inverse operator 

{L{v)) - \)-^ : L\R)^D{L{v)) 

is a continuous linear operator. Thus, 

(2.53) A(Q, A) [-2, 2] {L{v) - A)"^ : L^{W) D{L{v)) exists. 

Next we apply Lemma[3]with A = L[v) and B = L{v). Choose A G C so that 
A(Q,A) [-2,2] (in view of (|2.45p . see the remark before Theorem [U A{Q,X) 
is a non-constant entire function, so such a choice is possible). Then, in view 
of (j2.43p . we have that A{Q,X) [—2,2] also. In view of the above discussion, 
this means that the operator L{v) — A maps bijectively D{L{y)) onto L^(R) and 
L{y) — A maps bijectively D{L{v)) onto L^(M). Thus, by Lemma [3l D{L{v)) is 
dense in L^(M) and L{v)* = L{v), i.e., (a) and (b) hold. 

Finally, in view of ()2.53p . (c) follows readily from part (b) of Lemma [2j 

□ 

3. Theorem S] shows that the spectrum of the operator L{v) is described by 
the equation (j2.50p . As we are going to explain below, this fact implies that the 
spectrum Sp{L{v)) could be described in terms of the spectra of the operators 
Lq = Lq{v), 9 £ [0,7r], that arise from the same differential expression £ = £q 
when it is considered on the interval [0, vr] with the following boundary conditions: 

(2.54) y{n) = e''y{0), {y' - Qy)[7,) = e'\y' - Qy)(0). 
The domains D{Lg) of the operators Lq are given by 

(2.55) D{Le) = {y£H^ : y' - Qy & Wli[0,7r]), holds, i{y) £ H^} , 
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where 

H' = H\[0,n]), = L\[0,7t]). 

We set 

(2.56) Le{y) = i{y), y £ D{Le). 

Notice that if y G i?^([0,7r]), then (giy) = f G ^^([0,7r]) if and only if n = 
y' " Qy £ ^'^d the pair (y, n) is a solution of the non-homogeneous 

system ([2391) . 

Lemma 5. Let ( ) anc? ( ) he the solutions of the homogeneous system 



^Ui J \U2 

112. 20\) which satisfy 

(SDK? 

// 

(2.58) A = yi(7r) + U2(7r) / 2cos(9, 6'G[0,7r], 

then the non-homogeneous system i2.19\) has, for each f G H^, a unique solution 
{y,u) = (iii(/),i?2(/)) such that 

Moreover, Ri is a linear continuous operator from into , and R2 is a linear 
continuous operator in with a range in Wi{[0,Tr]). 

Proof. By the variation of parameters method, every solution of the non-homogeneous 
system (I2.19P has the form 

(!;!:.0-'"w(.!:;g)+-w(:g 

where 

/•X j-X 

(2.61) vi{x) = - y2{x)f{t)dt + Cu V2{x) = yi{x)f{t)dt + C2. 

Jo Jo 

We set for convenience 

(2.62) miif) = - r y2{t)f{t)dt, m2{f) = T yi{t)f{t)dt. 

Jo Jo 

By (l2:60]) - ([2:62]) . the condition ([239]) is equivalent to 

(2.63) + C) (^;«) + (,„.(/) + C.) (»«) = (g 

This is a system of two linear equations in two unknowns Ci and C2. The corre- 
sponding determinant is equal to 

det ( y^^^)-f' ^=l + e'^' - A . e^^ = e^^(2cos0 - A). 
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Therefore, if (j2.58p holds, then the system ()2.63p has a unique solution 



where Ci = Ci(/) and C2 = C2{f) are linear combinations of mi{f) and m2{f). 
With these values of Ci(/) and C2{f) we set 

^l(/) =Vl-yi+V2- 2/2, R2{f) =Vi-Ui+V2- U2. 

By (]2.6ip and (j2.62p . the Cauchy inequality implies 

\viix)\ < r \y2it)f{t)\dt + |C7i(/)| < A . 11/11, \v2{x)\ < B ■ ||/||, 
Jo 

where A and B are constants. From here it follows that Ri and R2 are continuous 
linear operators in . Since 

— R^{f)=Vi—+V2 — , R2if)=Vl—+V2— + f, 

ax ax ax ax ax 

it follows also that Ri acts continuously from into H^, and R2 has range in 
Ty/([0, vr]), which completes the proof. 

□ 

Theorem 6. Suppose v £ Hi^^{M)) is -k -periodic. Then, 

(a) for each 9 G [0, vr], the domain D{Lq{v)) G (|i^.55p is dense in H^; 
(h) the operator Lg{v) G ^2. 5(j^ is closed, and its conjugate operator is 

(2.64) Le{vy = Le{v). 

In particular, if v is real-valued, then the operator Lg{v) is self-adjoint. 

(c) the spectrum Sp{Lg{v)) of the operator Lg{v) is discrete, and moreover, 

(2.65) Sp{Lg{v)) = {X £ C : A(A) = 2cose}. 

Proof. Firstly, we show that the operators Lg{v) and Lgijj) are formally adjoint, 
i.e., 

(2.66) {Le{v)y, h) = (/, Le{v)h) if y G D{Lg{v)), h G D{Le{v)). 
Indeed, in view of (j2.54p . we have 

{Le{v)y,h) = - lQ{y)hdx = - I {-{y - Qy)' h - Qy h) dx 
Jo Jo 

= (y - Qy)h „ + - / (y - Qy)h'dx - / Qy'hdx 

= + - / (y'W - Qyh! - Qy'h) dx. 
Jo 



[y'h' - Qyh' - Qy'h) dx = [y, Lg{v)h) , 



The same argument shows that 

1 r 

TT Jo 

which completes the proof of (j2.66p . 

Now we apply Lemma [3] with A = Lg{v) and B = Lg{v). Choose A G C so 
that A(Q,X) 7^ 2cos0 (as one can easily see from the remark before Theorem U 
A{Q, A) is a non-constant entire function, so such a choice is possible). Then, in 
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view of ()2.43p . we have that A(Q, A) / 2 cos 6 also. By LemmaEl Lg{v) — A maps 
bijectively D{Lg{v)) onto and Lq{v) — A maps bijectively D{Lg{v)) onto H^. 
Thus, by LemmaEl D{L0{v) is dense in and Lg{v)* = Lg{v), i.e., (a) and (b) 
hold. 

If A(Q, A) = 2cos0, then e*^ is a root of the characteristic equation (|2.22|) . so 
there is a special solution {ip, -0) of the homogeneous system (I2.20p (considered 
with C = -A) such that holds with p = e'^. But then (p £ D{Lg{v)) and 

Lq{v)ip = Xif, i.e., A is an eigenvalue of Lg{v). In view of Lemma El this means 
that ()2.65p holds. Since A(Q,A) is a non-constant entire function (as one can 
easily see from the remark before Theorem \TT\i the set on the right in ()2.65p is 
discrete. This completes the proof of (c). □ 

Corollary 7. In view of Theorem and Theorem we have 

(2.67) SpiL{v))= U SpiLeiv)). 

ee[o,7r] 

In the self-adjoint case (i.e., when v, and therefore, Q are real-valued) the 
spectrum SpiLiv)) C M has a band-gap structure. This is a well-known result 
in the regular case where v is an Lf^^(M)-function. Its generalization in the 
singular case was proved by R. Hryniv and Ya. Mykytiuk 8]. 

In order to formulate that result more precisely, let us consider the following 
boundary conditions (be): 

(a*) periodic Per+ : ^(vr) = y(0), (y' - Qy) (vr) = (y' - Qy) (0); 

(b*) antiperiodic Per' : y(7r) = -y(0), (y' - Qy) (tt) = - (y' - Qy) (0); 

Of course, in the case where Q is a continuous function, Per^ and Per^ 
coincide, respectively, with the classical periodic boundary condition yin) = 
y(0), y'in) = y'(0) or anti-periodic boundary condition ^(vr) = — y(0), y'in) = 
—y'iO) (see the related discussion in Section 6.2). 

The boundary conditions Per^ are particular cases of (j2.59p . considered, re- 
spectively, for 9 = OT 6 = TT. Therefore, by Theorem [6l for each of these two 
boundary conditions, the differential expression (j2.18p gives a rise of a closed 
(self adjoint for real v) operator Lpf,^± in = L^([0,7r]), respectively, with a 
domain 

(2.68) DiLp,,+ ) = {yeH^ : y' - Qy G <([0, ^]), (a*) /loMs, /(y) G 
or 

(2.69) DiLp,,-) = {yeH^ : y' - Qy &Wli[^M), ih*) holds, /(y) G 

The spectra of the operators Lpg^± are discrete. Let us enlist their eigenvalues 
in increasing order, by using even indices for the eigenvalues of Lp^^+ and odd 
indices for the eigenvalues of Lp^j.- (the convenience of such enumeration will be 
clear later): 

(2.70) Sp (Lpe.r+) = {Ao, A2 , A^, A4 , A4 , Ag , Ag , . . .}, 

(2.71) 5p(Lpe,-) = {Ar,A+,A3,A+,A5,A+...}. 
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Proposition 8. Suppose v = C + Q' , where Q G Lf^^CR)) is a ir-periodic real 
valued function. Then, in the above notations, we have 

(2.72) Ao < A;f < A+ < A2 < < A3 < A+ < A4 < A| < Ag < < • • • . 

Moreover, the spectrum of the operator L{v) is absolutely continuous and has a 
band-gap structure: it is a union of closed intervals separated by spectral gaps 

(-cx),Ao), (A^,A^), (A2,A^),--- ,(A;;,A+),--- . 

Let us mention that A. Savchuk and A. Shkalikov [25| have studied the Sturm- 
Liouvihe operators that arise when the differential expression £q, Q S L^([0, 1]), 
is considered with appropriate regular boundary conditions (see Theorems 1.5 
and 1.6 in [27]^. 

3. Fourier representation of the operators Lp^^± 

Let L^^ denote the free operator = —(f'/dx^ considered with boundary 
conditions be as a self-adjoint operator in L^([0,7r]). It is easy to describe the 
spectra and eigenfunctions of L^^ for be = Per^, Dir : 

(a) Sp{L^p^^+) = {n^, n = 0,2,4, . . .}; its eigenspaces are ii^JJ = S'panje^*"^} 
for n > and Eq = {const}, dim.E'^ = 2 for n > 0, and dimi?Q = 1. 

(b) Sp{LPp^^_) = {n^, 71 = 1, 3, 5, . . .}; its eigenspaces are = 5pan{e^*"^'}, 
and dimE"^ = 2. 

(c) Sp{L^^^^) = {n^, n G N}; each eigenvalue is simple; a corresponding 
normalized eigenfunction is \/2sinnx. 

Depending on the boundary conditions, we consider as our canonical orthog- 
onal normalized basis (o.n.b.) in L^([0,7r]) the system Uk{x), k G r^c, where 

(3.1) if be = Per^ = ex.p{ikx), k € ^per+ = 2Z; 

(3.2) if be = Per~ Uk = exp{ikx), k £ Fpg^- = 1 + 2Z; 

(3.3) if be = Dir Uk = V^sinkx, k G F^jr = N. 

Let us notice that {uk{x), k G Fbd is a complete system of unit eigenvectors of 
the operator L^^. 

We set 
(3.4) 

H},^^^ = {f € : /(vr) = /(0)}, H},^^. = {f € : /(vr) = -/(O)} 
and 

(3.5) Hh,r = {/ G : /(^) = /(O) = O} . 

One can easily see that {e*^^, k G 2Z} is an orthogonal basis in Hpgj.+ , {e*^^, k G 
1-I-2Z} is an orthogonal basis in Hp^^_, and {^/2sinkx, k G N} is an orthogonal 
basis in Hp^^. 

From here it follows that 

(3.6) Hi = i fix) = ^ fkUkix) : = (1 + k')\fk\' < 00 i . 

I fcePic fcePtc J 
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The following statement is well known. 

Lemma 9. (a) If f,g G L\[0,it]) and f ~ Efce2Z /fce^''^' 9 ~ Efce2z5fce*''^ 
are their Fourier series respect to the system {e^^^, k G 2Z}, then the following 
conditions are equivalent: 

(i) f is absolutely continuous, /(vr) = /(O) and f'{x) = g{x) a.e.; 

(a) gk = ikfk VA; e 2Z. 



(b) Iff,g G Li([0,7r]) and f ~ Efc6i+2zAe^'"> 9 ^ Y.kei+2z9ke'^'' 



are 



their Fourier series respect to the system {e^^^ ,k G 1 + 2Z}, then the following 
conditions are equivalent: 

(i*) / is absolutely continuous, f[-K) = — /(O) and f'{x) = g{x) a.e.; 
{a*) gk = ikfk V/fcGl + 2Z. 

Proof. An integration by parts gives the implication (i) =^=- (ii) [or [i*) ^ 

To prove that (ii) =^ (i) we set G{x) = g{t)dt. By (ii) for /c = 0, we have 
G(7r) = g{t)dt = irgo = 0. Therefore, integrating by parts we get 

gf^ = - r g{x)e-'^''dx = - r e-'''''dG{x) = ikGk, 
Jo Jo 

where Gk = ^ Jq e~^^^G{x)dx is the A;-th Fourier coefficient of G. Thus, by (ii), 
we have Gk = fk for k ^ 0, so by the Uniqueness Theorem for Fourier series 
f{x) = G{x) + const, i.e.(i) holds. 

Finally, the proof of the implication (ii*) ^ (i*) could be reduced to part 
(a) by considering the functions f{x) = f{x)e'^^ ~ Z]fcei+2Z /fc-i^*'^^ ^"^^ 9i^) ~ 
g{x)e^^ + z/(x)e*^. We omit the details. □ 

The next proposition gives the Fourier representations of the operators Lp^^± 
and their domains. 

Proposition 10. In the above notations, if y ^ Hp^^±, then we have y = 
Y^r + Vk^'^^ ^ D{Lpf,r±) and t{y) = h = Xlr + hkC^^^ G if and only 
if 

(3.7) hk = hk{y) := k'^yk+ ^ V{k - m.)ym + Cyk, ^|/ifc|^<oo, 



I.e., 



(3.8) D{Lp,,±) = {y G F],^,± : {hk{y))ker^^^± e f (Fp^.i)} 
and 

(3.9) Lp,My)= E 

Proof. Since the proof is the same in the periodic and anti-periodic cases, we con- 
sider only the case of periodic boundary conditions. By (|2.68p . if y G D{Lp^^+), 
then y G H\,^^+ and 

(3.10) l{y) = -z' -Qy' + Cy = he L\% vr]), 
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ikx 



where 

(3.11) z:=y' -QyeW^{[0,7r]), ^(vr) = z(0). 
Let 

ke2Z fee2Z fce2Z 

be the Fourier series of y,z and h. Since ^(vr) = z{0), Lemma [9] says that the 
Fourier series of z' may be obtained by differentiating term by term the Fourier 
series of z, and the same property is shared by y as a function in Hp^^+. Thus, 
([Hl^ imphes 

(3.12) - ikzk - ^ - m)imym + Cyk = 

m 

On the other hand, by p. lip , we have z^ = iky^ — q{k — m)ym, so substi- 
tuting that in (j3.12p we get 



(3.13) ~ik 



m)yr, 



qik - m)imym + Cyk = hk, 



which leads to (|3.7p because V{m) = imq{m), m £ 2Z. 

Conversely, if ()3.7p holds, then we have ()3.13p . Therefore, (j3.12p holds with 
Zk = ikyk - J2m 'ii^ - '^)ym- 

Since y = YlVk^^^^ G Hp^^+, the Fourier coefficients of its derivative are 
ikyk, k G 2Z. Thus, (z^) is the sequence of Fourier coefficients of the function 
z = y' -QyGLH[0,7r]). 

On the other hand, by (I3.12p . (ikz^) is the sequence of Fourier coefficients 
of an L^([0, 7r])-function. Therefore, by Lemma [U the function z is absolutely 
continuous, ^(vr) = -2(0), and {ikzk) is the sequence of Fourier coefficients of its 
derivative z' . Thus, (|3.10p and (j3.1ip hold, i.e., y G D{Lpf,^+) and Lp^^^+y = 
l{y) = h. □ 

Now, we are ready to explain the Fourier method for studying the spectra of 
the operators Lper±- L^t 

be the Fourier isomorphisms defined by corresponding to each function / G 
the sequence (/fc) of its Fourier coefficients fk = (/, "Ufe), where {u^, k G Fpg^i} is, 
respectively, the basis dsn) or ([32]). Let T'^ be the inverse Fourier isomorphism. 
Consider the unbounded operators >C+ and C- acting, respectively, in £^(rpg,,±) 

as 

(3.14) C±{z) = {hk{z))^^^ , hk{z) = k^Zk+ Y] V{k - m)zm + Czk, 
respectively, in the domains 

(3.15) D{L^) = {z^e{\k\,Tp^,±): £±(z) G ^^(rpg.i)} , 
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where ^^(|A;|, Fpg^-i) is the weighted £^-space 



In view of (|3.6p and Proposition 1101 the following theorem holds. 
Theorem 11. In the above notations, we have 

(3.16) D{Lp,,±) = J^-^ {D{C±)) 
and 

(3.17) Lp^r± = T'^ oC±oT. 

If it does not lead to confusion, for convenience we will loosely use one and 
the same notation Lp^^± for the operators Lp^^± and C±. 

4. Fourier representation for the Hill-Schrodinger operator with 

dirichlet boundary conditions 

In this section we study the Hill-Schrodinger operator L£)ir{v), v = C + Q' , 
generated by the differential expression ^Q(y) = —{y' — Qy)' — Qy' considered on 
the interval [0, vr] with Dirichlet boundary conditions 

Dir : y{0) = ^(vr) = 0. 

Its domain is 
(4.1) 

D{LdUv)) = {y e : y' -QyeWl{[0,Tr]), y(0) = ^(vr) = 0, iqiy) e H^^} 
and we set 

(4.2) LDir{v)y = £Q{y). 

Lemma 12. Let j '^''^d i^j solutions of the homogeneous system 



\2. 20\} which satisfy 
(4.3) 



yi(o)\ ^ (i\ (y2io)\ ^ (0 

ni(0)j [oh U(0)j ll 



// 

(4.4) y2(vr) / 0, 

then the non-homogeneous system i2.19\) has, for each f G H^, a unique solution 
{y,u) = (i?i(/),ii2(/)) such that 

(4.5) y(0) = 0, y{7T) = 0. 

Moreover, Ri is a linear continuous operator from into H^, and R2 is a linear 
continuous operator in with a range in WI{\}},'k\). 
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Proof. By the variation of parameters method, every solution of the non-homogeneous 
system (I2.19P has the form 



^u[x) J yui[x) J yu2(x) J 

where 

(4.6) vi{x) = - y2{x)f{t)dt + Ci, V2{x)= yi{x) f {t)dt + C2. 



By (I4.3p . the condition y(0) = wih be satisfied if and only if Ci = 0. If so, the 
second condition y(7r) = in (14. Sp is equivalent to 

mi(/)yi(vr) + (m2(/) + C2)y2(vr) = 0, 

where 

^lU) = - r y2{x)f{t)dt, m2{f) = r yi{x)f{t)dt. 



Thus, if y2(7r) 7^ 0, then we have unique solution (y,n) of (12.190 that satisfies 
(14. 5p . and it is given by (14. 6p with Ci = and 



(4.7) C2(/) = -^mi(/)-m2(/). 

2/2 (tt) 

Thus, we have I ^I^n I = ( {^"'^I'^l I , where 
\u{x)} \R2{f)l 



RiU) ={ - I y2{x)f{t)dt ] ■ yiix) + ( / yi{x)f{t)dt + C2{f) ] ■ y2{x) 

and 

R2{f) ={ - I y2{x)f{t)dt ) • ni(x) + ( / yi{x)f{t)dt + C2{f) ) • U2{x). 



It is easy to see (compare with the proof of Lemma [S]) that Ri is a linear con- 
tinuous operator from into H^, and R2 is a linear continuous operator in 
with a range in Wl ( [0, vr] ) . We omit the details. □ 



Now, let us consider the systems (j2.19p and (|2.20p with a spectral parameter 

(v^{x X)\ f V'){x X)\ 

, \{\ and ( / ' x\ I be the solutions 
ui{x,\)) \U2{X,X)J 

of the homogeneous system (|2.2Up that satisfy (|4.3p for x = 0. Notice that 



(4.8) y2iv;x,X) = y2{v;x,X) 

hoc 



Theorem 13. Suppose v G if^^^(M) is -k -periodic. Then, 



(a) the domain D{L£)ir{v)) S d^. jp is dense in H^; 
(h) the operator L^nr^v) is closed, and its conjugate operator is 

(4.9) (LDiAv))* = LDiriv). 

In particular, if v is real-valued, then the operator Loiriy) is self-adjoint. 
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(c) the spectrum Sp{LDir{v)) of the operator L£)ir{v) is discrete, and moreover, 

(4.10) Sp{LDir{v)) = {A G C : y2(vr, A) = 0}. 

Proof. Firstly, we show that the operators L£,ir{v) and Ljjiriv) are formally 
adjoint, i.e., 

(4.11) {LDir{v)y,h) = {f,LDir{v)h) ii y e D{LDir{v)), he D{LDir{v)). 
Indeed, in view of ()4.ip . we have 

{LDir{v)y, h) = - I £Q{y)hdx = - {-{y - Qy)'h - Qy'h) dx 
^ JO JO 

= — (y - Qy)h n + - / (y - Qy)h'dx - / Qy'Mx 

71" " TT 7o Jo 

= + - I iy'hl -Qyhl -Qy'h)dx. 
Jo 



The same argument shows that 

1 r 

71" io 



y'h' - Qyh' - Qy'h) dx = (y, LDir{v)h) 



which completes the proof of (|4.11|) . 

Now we apply Lemma [3] with A = LDir{v) and B = LDir{v). Choose A G C 
so that y2{v;7r,X) / (in view of ()2.44p . see the remark before Theorem [TTl 
y2{v]'K,\) is a non-constant entire function, so such a choice is possible). Then, 
in view of (|4.8p . we have y2{y;T^,X) ^ also. By Lemma [121 LDir{v) — A maps 
bijectively D{L£iir{v)) onto and L£)ir{v) — X maps bijectively D{L£)ir{v)) onto 
H^. Thus, by Lemma [3l D{LDir{v)) is dense in and {L£)ir{v))* = Lfjiriv), 
i.e., (a) and (b) hold. If y2{v;TT,X) = 0, then A is an eigenvalue of the operator 
LDiriv), and y2iv; x, A) is a corresponding eigenvector. In view of Lemma [T2| this 
means that that (j4.10p holds. Since y2{'K,\) is a non-constant entire function, 
the set on the right in (I4.10p is discrete. This completes the proof of (c). □ 

Lemma 14. (a) If f,g G L^{[0,tt]) and f ~ ^^^^ /fe\/2 sin fcx, g ^ go + 
Yl'k^=i SkV^ cos kx are, respectively, their sine and cosine Fourier series, then 
the following conditions are equivalent: 

(i) f is absolutely continuous, /(O) = /(tt) = and g{x) = f'{x) a.e.; 
(a) 50 = 0, gk = kfk V/c G N. 

(b) Iff,g G Li([0,7r]) andf ~ fo+Y.kLi fkV^ cos kx and g ~ Y.kLi 9kV2 sin kx 
are, respectively, their cosine and sine Fourier series, then the following condi- 
tions are equivalent: 

(i*) f is absolutely continuous and g{x) = f'{x) a.e.; 
(n*) 9k = -kfk k£N. 

Proof (a) We have (i) {ii) because go = I fo 9{x)dx = ;^(/(vr) - /(O)) = 0, 
and 

gk = — g{x)V2coskxdx = —f{x)\/2coskx H — / f{x)V2s\nkxdx = kfk 
Jo vr Jo 
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for every A; G N. 

To prove that (ii) ^ (i), we set G{x) = /q g{t)dt; then G{tt) = G{0) = 
because go = 0. The same computation as above shows that gk = kGk V/c G N, 
so the sine Fourier coefficients of two L^-functions G and / coincide. Thus, 
G{x) = f{x), which completes the proof of (a). 

The proof of (b) is omitted because it is similar to the proof of (a). □ 

Let 

oo 

(4.12) Q^Yl q{k)V2sm kx 

k=l 

be the sine Fourier expansion of Q. We set also 

(4.13) V{0) = 0, V{k) = kq{k) for A; G N. 

Proposition 15. In the above notations, ify G Hjjir, then we have y = X^fcLi Vk sin /ex G 

D{L£)ir) and £{y) = h = X^fcLi hkV^sinkx G if and only if 

(4.14) 

^ oo 

hk = hk{y) = k^yk+^ ^ {y{\k - m\) - V{k + m)^ ym+Cyk, ^ |/ifcp < oo, 

^ k=l 



I.e., 



(4.15) 

oo 

D{LDir) = {y ^ H},^, : (/ifc(y)r G^'(N)}, Lz)i.(y) = ^ /ifc(y) ^2 sin fca;. 

k=l 

Proof. By gl]), if y G D{LDir), then y G Hjy^^ and 

i{y) = -z' -Qy' + Cy = h€L^{[0,7T]), 

where 

(4.16) z:=y' -QyeWl{[0,7r]). 

Let 

oo oo oo 

y ^ ''^^ykV2smkx, z ~ z^V^ cos /cx, /i ~ /ifc\/2 sin /cx 

k=l k=l k=l 

be the sine series of y and /i, and the cosine series of z. Lemma [HI yields 

oo oo 

z' ~ ^^(— /czfc)\/2 sin /cx, y ^ kykV2cos kx. 

k=l k=l 

Therefore, 



(4.17) hk = kzk - {Qy')k + Cyk, A; G N, 

where {Qy')k are the sine coefficients of the function Qy' G L-'^([0, vr]). 
By (jCrej) . we have 

Zk = kyk - {Qy)k, 
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where {Qy)k is the k-th cosme coefficient of Qy. It can be found by the formula 

iQy)k = - Q{x)y{x)V2 cos kxdx = • y^, 

71" Jo 1 

m=l 



with 



1 /"'^ 1 

(k) = — / (5(x)\/2 cos A;x\/2 sinmxdx = — / Q(x)[sm{m+k)x+sm{m—k)x]dj 



I: 



g(m + k) + q{m — k), m > k 
^ X q{2k), m = k 

\q{m + k) — q{k — m) m < k. 



Therefore, 
(4.18) 



fe-i 



(Qy)fc = ^ X] ^("^ + ~ 75 X] ~ + 7f X] ^("^ ~ ^) 

^ m=l ^ m=l ^ 



m=A;+l 

In an analogous way we can find the sine coefficients of Qy' by the formula 



{Qy')k = — Q{x)y' {x)V2s'mkxdx = h 
^■^0 ^1 



m ' '^y-mi 



where hm are the cosine coefficients of (5(x)\/2 sin /cx, i.e., 



bm{k) 



Q{x)V2sinkxV2 cos 



mx 



TT 



Q{x) [sin(A;+m)3;+sin(A:— m)x](ix 



1 

71 



q{k + m) + g(/c — m), m < k, 
q{2k), m = k, 

q{k + m) — g(m — fe) m > k. 



Thus we get 
(4.19) 



fc-i 



{Qy')k = -^"^ q{m+k)mym + ^'^ q{k-m)mym--^ ^ q{m-k) 

^ m=l ^ m=l ^ 



myv 



Finally, KW\ and (f4T9]) . imply that 



^^y/fc - ^(Qy)fc - (Qy') 



m,=fc+l 



fc-1 



k^yk- 



(rn-+fc)g(m+A:)-| — p {m—k)q{m—k)-\ — p ^^^{k—ni)q{k- 

m=l m=fc+l m=l 

Hence, in view of (|4.13p . we have 

/ifc = k'^yk + y^Yl (^(1^ ~ "^D ~ '^(^ + "^)) 2^"^ + '^y^' 

^ )n=l 

i.e., KIM holds. 
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Conversely, if ()4.14p holds, then going back we can see, by ()4.17p . that z = 
y' ~ Qy £ -^^([0) ^]) has the property that kzk, A; e N, are the sine coefficients of 
an L^([0, 7r])-function. Therefore, by Lemma \T4\ z is absolutely continuous and 
those numbers are the sine coefficients of its derivative z' . Hence, z = y' — Qy S 
l^/([0,7r]) and £{y) = h, i.e., y G D{LDir) and LDiriv) = h. □ 

Let 

T-.H"^^ e{N) 

be the Fourier isomorphisms that corresponds to each function / € the se- 
quence {fk)k&^ of its Fourier coefficients fk = (/, \/2sin/cx), and let J-~^ be the 
inverse Fourier isomorphism. 

Consider the unbounded operator and acting in ^^(N) as 
(4.20) 

Cd{z) = (hkiz)))^^^ , hk{z) = k'^Zk+—i= ^ {y{\k - m\) - V{k + m)^ z^+Czk 
in the domain 

(4.21) D{Ld) = {z<-e{\k\,^)- Ld{z)^e{n)], 
where ^^(|/c|,N) is the weighted ^^-space 

^'(|A:|,N) = |z = (zfc)fceN : ^ \A\k? < ooj . 

In view of (|3.6p and Proposition 1151 the following theorem holds. 
Theorem 16. In the above notations, we have 

(4.22) D{LDir) = {D{Cd)) 
and 

(4.23) Loir = oCdoT. 

If it does not lead to confusion, for convenience we will loosely use one and 
the same notation Lcur for the operators LDfr and Cd- 

5. Localization of spectra 
Throughout this section we need the following lemmas. 
Lemma 17. For each n G N 



(5.1) < 



1 2 log 6n 
< 2 — • 

— fe^l n ' 



^14 

The proof is elementary, and therefore, we omit it. 
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Lemma 18. There exists an absolute constant C > such that 
(a) if n G N and 6 > 2, then 



^ ,n--/c2| + 5 - ^ ' 

(b) if n >0 and b > then 

1 C 

|„2_p|2 + 52 ^ („2 + 52) 1/2(^4 + 52) 1/4- 

A proof of this lemma can be found in jjJj, see Appendix, Lemma 79. 

We study the locahzation of spectra of the operators Lp^^± and L^ir by using 
their Fourier representations. By (I3.14p and Theorem II 1^ each of the operators 
L = Lpf,j.± has the form 

(5.5) L = L^ + V, 

where the operators and V are defined by their action on the sequence of 
Fourier coefficients of any y = X^r^ ± Vk 6xp ikx S Hp^^^± : 

(5.6) L": {y^)^{ey^), keTp,,± 
and 

(5.7) V: {ym)^izk), z,, = J^F(A; 

m 

(We suppress in the notations of and V the dependence on the boundary 
conditions Per^.) 

In the case of Dirichlet boundary condition, by ()4.20p and Theorem 1 161 the op- 
erator L = L has the form (j5.5p , where the operators and V are defined by 
their action on the sequence of Fourier coefficients of any y = ykV^sinkx G 

(5.8) lO : (yfc) ^ (k^yk), keN 
and 

(5.9) ^ 

^:(ym)^(^fe), Zfe = — =^ (yd/c-ml) -y(/c + m)) A;,mGN. 

(We suppress in the notations of and V the dependence on the boundary 
conditions Dir.) 

Of course, in the regular case where v G ^^([0, vr]), the operators and V 
are, respectively, the Fourier representations of —d^/dx^ and the multiplication 
operator y ^ v ■ y. But if ?; G H^^{M.) is a singular periodic potential, then 
the situation is more complicated, so we are able to write (15. 5p with (15.61) and 
(15. 7p . or (15. Sp and (|5.9p . only after having the results from Section 3 and 4 (see 
Theorem [TT] and Theorem I16p . 
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In view of ()5.6I) and (15. 8p the operator LP is diagonal, so, for A / /c^. A; G Ti,c, 
we may consider (in the space I'^iTi)^)) its inverse operator 

(5.10) Rl : (zfe) ^ (y^) , ^ e U,. 

One of the technical difficulties that arises for singular potentials is connected 
with the standard perturbation type formulae for the resolvent Rx = {\ — LP — 
V)^^ . In the case where v E -Z>^([0, tt]) one can represent the resolvent in the form 
(e.g., see [1], Section 1.2) 



oo 

(5.11) Rx = {l- Rlvr'Rl = Y,{RlvfRl 

k=0 



or 

oo 

(5.12) Rx = R'xil - VRlr' = Rl{VRl)\ 

k=0 

The simplest conditions that guarantee the convergence of the series (j5.1ip or 
(f5l2]) in f are 

\\RxV\\ < 1, respectively, \\VRI\\ < 1. 

Each of these conditions can be easily verified for large enough n if Re A G 
[n — l,n + 1] and |A — > C(||f ||), which leads to a series of results on the 
spectra, zones of instability and spectral decompositions. 

The situation is more complicated if t; is a singular potential. Then, in general, 
there are no good estimates for the norms of Rx^ and VR^. However, one can 
write ([5TT]) or ([5T2]) as 

oo 

(5.13) Rx = Rl + R^VRl + RlVR^VRl + ... = kI+Y, KxiKxVKxTKx, 

m=l 

provided 

(5.14) {Kx? = Rl 

We define an operator K = Kx with the property (j5.14p by its matrix represen- 
tation 

(5.15) Kjm = p^TT^^yiT^'^im' i'"^ S Tbc, 
where 

= ^e*'^/2 z = re''^, < (/? < 27r. 
Then Rx is well-defined if 

(5.16) WKxVKx : t\Tkc) ^ ^'(r^c)!! < 1- 

In view of (I2.14p . ()5.7p and (I5.15p . the matrix representation of KVK for 
periodic or anti-periodic boundary conditions he = Per^ is 

V{j — m) i{j — m)q{j — m) 



(5.17) (KVK) 



(A - j2) V2(A - m2) 1/2 (A - j2)l/2(A - ^2)1/2 
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where j, m E 2Z for be = Per^, and j,m £ 1 + 2Z for be = Per~ . Therefore, we 
have for its Hilbert-Schmidt norm (which majorizes its ^■^-norm) 



(5.18) WKVKWls = E 



By (I4.13p . (j5.9p and (|5.15p . the matrix representation of KVK for Dirichlet 
boundary conditions be = Dir is 



(5.19) (KVK) 



1 V{\j-m\) 1 V{j + m) 



V2(A- j2)V2(A-m2)V2 ^(A-j2)i/2(A_m2)i/2 
\j - m\q{\j - m\) 1 (j + m)g(j + m) 



V2 (A - j2)i/2(A - m2)i/2 ^ (A - j2)i/2(A - m2)i/2 ' 

where j, m G N. Therefore, we have for its Hilbert-Schmidt norm (which ma- 
jorizes its ^^-norm) 
(5.20) 

~ ^ |A - i2||A - m2| ^ |A - 72||A - m2| 

j,mGN ' " ' j,meN ' " ' 

We set for convenience 
(5.21) g(0)=0, r{s)=q{\s\) for s / 0, s € Z. 

In view of (j5.20p and (j5.2ip . we have 

(j — m)2|f(j — m)p 



(5.22) \\KVK\\ls< 



A — ?'2| |A — m^l 

We divide the plane C into strips, correspondingly to the boundary conditions, 
as follows: 

if be = Per+ then C = HqU H2U H4U ■ ■ ■ , and 
if be = Per- then C = HiU H3U H5U ■ ■ ■ , 
where 

(5.23) i/o = {A e C : iJeA < 1}, Hi = {X £ C : Re X < A}, 

(5.24) Hn = {XeC : {n-lf < ReX< {n + lf}, n > 2; 
- if 6c = Dir, then C = Gi U G2 U G3 U • • • , where 

(5.25) Gi = {X: ReX< 2}, G„ = {X : {n - l)n < Re X < n{n + 1)}, n>2. 
Consider also the discs 

(5.26) L>n = {A e C : |A - n^l < n/4}, n G N, 
Then, for n > 3, 

fcen+2Z ' ' fcen+2Z ' ' 
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and 

fcez ' ' fcez ' ' 

where Ci is an absolute constant. 

Indeed, if A G Hn, then one can easily see that 

|A - fc^l > In^ - /cV4 for ken + 2Z. 
Therefore, if A G Hn \ Dn, then ()5.ip imphes that 

y 1 < ^ I y 4 ^ 8 ^ 81og6n ^^^logn 

which proves the first inequality in (|5.27|) . The second inequality in (|5.27|) and 
the inequalities in ()5.28p follow from Lemma [T7] by the same argument. 

Next we estimate the Hilbert-Schmidt norm of the operator KyyK\ for he = 
Per^ or Dir, and correspondingly, A G Hn \ Dn or A G G„ \ Dn, n G N. 

For each £^-sequence x = {x{j))j,z^ and m G N we set 

(5.29) £.^{x)= I \x{j)\'' 

\\j\>m 

Lemma 19. Let v = Q', where Q{x) = Ylke2z'ii^)^^^^ — Sm=i q{'m)^/2sinmx 
is a TT-periodic L^([0,7r]) funetion, and let 

q = {q{k))ke2Z, q = {q{m))mm 

be the sequences of its Fourier coefficients respect to the orthonormal bases {e^^^ , k G 
2Z} and {\/^smmx, m G N}. Then, for n > 3, 

(5.30) \\KxVKx\\HS<c(£^{q) + \\q\\/^y A G i7„ \ Z)„, 6c = Per±, 
and 

(5.31) \\KxVKx\\HS<c{£^{q) + \\q\\/V^^, A G G„, \ ^n, he = Dir, 
where C is an absolute constant. 

Proof. Fix n G N. We prove only (j5.30p because, in view of (j5.2ip and (j5.22p . the 
proof of (j5.3ip is practically the same (the only difference is that the summation 
indices will run in Z). 
By dSHI), 



s^ 



(5.32) lli.yi.lll,, < E E |A-^^||A-(n. + .)^| l^^^)'' = + + ^ 
where s G 2Z, m G n + 2Z and 

(5.33) Si= E ••• ' ^2= E ' ^3= E ■ 

\s\<y/n y/n<\s\<An |s|>4n 



3, 
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The Cauchy inequality implies that 

(5-34) ^ _ X _(m -L c^2| - X] 

men+2Z ' 

Thus, by (lOHl) and (lOOl) . 



|A — m^llA — (m + s)^| ' |A — m^p 

men+2Z ' " ^ ^ ' men+2Z ' ' 



(5.35) < i^f^hf = —M^ ^^Hn\D„ 

|s|<v^ 



(5.36) S2<(4n)2^ ^ |(7(s)|2 = 16Ci (£:^(g)) ' , A G \ 

Next we estimate S3 for n > 3. First we show that if |s| > An then 

(5.37) ^- f | <16^^, AGF„\I?.. 

|A — m^||A — (m + sj^l n 

Indeed, if |m| > |s|/2, then (since |s|/4 > n > 3) 

|A - m^l > - |i?e A| > s^/A -{n + if > - (|s|/4 + 1)^ > 

Thus, by dOT]) . 

8 ^g Cjlogre 

|A - m^llA - (m + s)2| ~ ^ |A - (m + s)2| ~ n 

|m|>|s|/2 

for A G i/„, \ Z)„. If \m\ < \s\/2, then |m + s| > |s| — \s\/2 = \s\/2, and therefore, 

|A - (m + s)2| > (m + sf - \Re A| > s'^/i - (n + 1)^ > s'^/8. 
Therefore, by ([OTl) . 

|A — m^llA — (m + s)-^| "■^^lA — m^l ~ n 

|m|<|s|/2 ' II V ' ' m ' ' 

for A G -ff„, \ Dn , which proves (|5.37p . 
Now, by (|537D . 

(5.38) S3 < 16^1^ |,(.)|^ = 16^1^ i£,U,)f . 

|s|>4n 

Finally, (|532D . (|535D . (|536D and (OH]) imply ([QUI) . 

Let denote the half-plane 

(5.39) H^ = {X£C: Re\< + N}, iV G N, 
and let Rn he the rectangle 

(5.40) i?7V = {AGC: -iV < i?e A < iV^ + iV, |/mA| < iV}. 



□ 
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Lemma 20. In the above notations, for be = Per^ or Dir, we have 

(5.41) sup{\\K^VK^\\hs, XeH^\RN}<c(^ ^^"^v^' 11^11 + ^a^Q)^ ^ 

where C is an absolute constant, and q is replaced by q ifbc = Dir. 

Proof. Consider the sequence r = {r{s))s<^z, defined by 

for odd s, 

max{\q{s)\,\q{—s)\) for even s, 



(5.42) r(s) 



Then, in view of (I5.42p . we have r € £^(Z) and ||r|| < 2\\q\\. 
If be = Per^, then we have, by (15.180 . 

(5.43) nKVKf„,< ^ 

On the other hand, if be = Dir, then (j5.22p gives the same estimate for ||iiryi<'|||^^ 
but with r replaced by the sequence f £ (|5.2ip . So, to prove (I5.4ip . it is enough 
to estimate the right side of (|5.43p for A G H'^ \Rn}. 
If ReX< -N, then ([S^IH]) imphes that 

\\kvk\\hs< 2^ ^vTjWT^' 

On the other hand, for 6 > 1, the following estimate holds: 

v-^ (?' — m)^|r(i — m)P ,, ,,r,l+7r 

5-44 y ^775 :^<^rf—^. 

|62+j2||52+^2| 5 

Indeed, the left-hand side of (j5.44p does not exceed 

2(j2 + m'^)\r{j — m)p 



E 



|62+j2||62 + 



^ 2 



m ' ' j j ' ' m 

< 4||r||^ ( i + 2 /" -TT^dx] = 4||r||2 ( -L + ^ ) < 4||r-||2ii:^. 



62 7o 6^ + a:2 7 ^ 
Now, with fe = ViV, (fOil) yields 



(5.45) \\KVKfHs<C^^ if ReX<-N, 



where C is an absolute constant. 

By (I5.43P and the elementary inequality 



|A — m^l = Y^(x — m2)2 + y2 > (^i^, _ rn'^\ -)_ |y|)/i/2j A = x + 
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we have 

(5.46) ||i^AFi^A||?/5<E^(^'y;")l^(") 

where 



(5.47) a{x,y;s) = Y, 



mez ^'^ - "^^1 + |y|)(k - + s)2| + |y|) ■ 
Now, suppose that X = x + iy £ \ Rn and \y\ > N. By (fOSj) and (|09]) . 

l<n<N 

SO A G for some n < N. Moreover, 

(5.48) a{x,y;s) <l6a{n'^,N;s) if A G G„, |y| > iV. 
Indeed, then one can easily see that 

|x - m^l + |y| > ^(|n^ - m^l + A^), m G Z, 

which imphes (j5.48p . 

By (15:171) and ([5A8\\ . ii X = x + iy £ Gn\ Rn and \y\ > N, then 

(5.49) ||i^Al^ifA|||^5< < + + 
where 

Si= J] cT(n2,iV;.)|r(s)|2, S2 = S3 = • • • . 

|s|<4y7V 4//V<|s|<4n |s|>4n 

If |s| < 4\/iV, then the Cauchy inequality and (|5.4[) imply that 

1 G '^2G 
a(r?, N; s) < 32N ■ V — ^ < 32N—- -— < 

^ ' ^ |n2 -m2|2 + Ar2 - Ar(n4 + iV2)i/4 - ^ 

Thus 

(5.50) Si<^||r||2. 



If 4viV < |s| < 4n then the Cauchy inequality and ()5.4p yield 

1 C 
cj(n^ N; s) < ?>2r? ■ V -^5 ^7^5 — r < 3271^—- — — < 32G 

^ ' ' ^- |„2 _^2|2 _^^2 - Ar(„4 _^ ^2)1/4 - 

because n < N. Thus 

(5.51) E2<32C-(^4^(r))'. 

Let |s| > 4n. If \m\ < \s\/2 then |m + s| > |s|/2, and therefore, 

In^ - (m + s)^| > |m + sj^ - > {\s\/2f - {\s\/Af > s'^/8. 



34 PLAMEN DJAKOV AND BORIS MITYAGIN 

Thus, by (IO]l . 

\m\^s\/2 N){\n'^ - (m + s)2| + iV) - ^ |n2 - m2| + iV - ^ ' 

If \m\ > \s\/2, then we have the same estimate because m? — > (|s|/2)2 — 
(|s|/4)2 > s'^/8, and therefore, again by (|5.3p . 

^ (|n2 - m2| + iV)(|n2 - (m + s)2| + iV) ~ ^ |n2 - (m + s)2| + iV ~ a/ZV ' 

|m|>|s|/2 ^' I /\i \ / I / ^ I \ y I V 



Thus C7(n2, AT; s) < 32C(log N)/VN, so we have 

(5.52) S3 < 32C||rfi^. 

Now, in view of (fOSI) and (fOT]) . the estimates (fOSj) and ([530]) - ([532]) yield 
(j5.4ip . which completes the proof. □ 

Theorem 21. For each periodic potential v G Hj^^{M), the spectrum of the 
operators Li,c{v) with be = Per^, Dir is discrete. Moreover, if be = Per^ then, 
respectively, for each large enough even number > or odd number N~ , we 
have 

(5.53) SpiLp,,±) cRn±U U 

n6Af±+2N 

where Rjy is the rectangle |5.^0[ ), Dn = {A : |A — n^j < n/4}, and 
(5.54) 

{2 /v+ -I- 1 
, # {Sp {Lp,r±) n Dn) = 2forne iV±+2N, 

where each eigenvalue is counted with its algebraic multiplicity. 
If be = Dir then, for each large enough number E N, we have 

00 

(5.55) Sp{LDn^)CRNlJ [j Dn 

n=N+l 

and 

(5.56) # {Sp {Loir) nRN)=N + l, # {Sp {Loir)) n D„) = 1 for n>N. 

Proof. In view of (j5.13p . the resolvent Rx is well defined if < 1. There- 

fore, (j5.53p and (j5.55p follow from Lemmas [19] and [20l 

To prove (j5.54p and (j5.56p we use a standard method of continuous parametriza- 
tion. Let us consider the one-parameter family of potentials v-r{x) = tv{x), t € 
[0, 1]. Then, in the notation of Lemma [191 have Vr = t ■ Q' , and the assertions 
of Lemmas 1191 and 1201 hold with q and q replaced, respectively, by r • g and r • q. 
Therefore, (I5.53P and (|5.55p hold, with L^c = Lbc{v) replaced by Lbc{vr). More- 
over, the corresponding resolvents Rx{Li,c{vt)) are analytic in A and continuous 
in r. 
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Now, let us prove the first formula in (I5.54p in the case be = Per^. Fix an 
even A^^ G N so that (j5.53p holds, and consider the projection 

(5.57) pN^^) = J_l {X-Lp,,^{vr))-'dX. 

^TTi Jx&dRN 

The dimension dim (P^(t)) gives the number of eigenvalues inside the rectangle 
Rn- Being an integer, it is a constant, so, by the relation (a) at the begging of 
Section 3, we have 

dimP^(l) = dimP^(O) = 2N+ + 1. 

In view of the relations (a)-(c) at the begging of Section 3, the same argument 
shows that (|5.54p and (15.56P hold in all cases. □ 

Remark. It is possible to choose the disks Dn = {A : |A — < r„} in 
Lemma[T9]so that r„/n 0. Indeed, if we take r„ = n/ip{n), where ip{n) — > oo 
but i~p{n)l^fn — > and if(n)E^{W) — > 0, then, modifying the proof of LemmafT^l 
one can get that ||i^A^-f^A||_f/5 — > as n — > oo. Therefore, Theorem [2T] could 
be sharpen: for large enough and N, i5.53\) - ^5. 56\) hold with Dn = {X : 
|A — n^l < r„} for some sequence {r„} such that rn/n 0. 

6. Conclusion 

The main goal of our paper was to bring into the framework of Fourier method 
the analysis of Hill-Schrodinger operators with periodic //^~^(R) potential, con- 
sidered with periodic, antiperiodic and Dirichlet boundary conditions. As soon as 
this is done we can apply the methodology developed in [151 El E] (see a detailed 
exposition in [4]) to study the relationship between smoothness of a potential 

V and rates of decay of spectral gaps 7n = A+ — A~ and deviations 5„ under 
a weak a priori assumption v G H~^. (In [151 El El |4j the basic assumption is 

V € i^^([0, vr]).) Still, there is a lot of technical problems; we present all the details 
elsewhere. But now let us give these results as stronger versions of Theorems 54 
and 67 in [4J. 

Theorem 22. Let L = + v{x) be a Hill-Schrodinger operator with a real- 
valued TT-periodic potential v G H^^(K), and let 7 = (7^) be its gap sequence. If 
uj = {uj{n))n£Z is a sub-multiplicative weight such that 

(R 1^ logu;(n) 

(6.1) \ as n —>■ 00, 

n 

then, with 

n = (n(n)), n(n) = '^^, 

n 

we have 

(6.2) -I ef{n,n)^v e H{n). 

If Q is a sub-multiplicative weight of exponential type, i.e., 

(6.3) hm > 0, 

n^oo n 
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then there exists £ > such that 

(6.4) 7 G£2(N,f)) ^ ?; G i7(e^l'"l). 

The following theorem summarizes our results about the Hill-Schrodinger op- 
erator with complex-valued potentials v G H~^. 

Theorem 23. Let L = LP + v{x) be the Hill-Schrodinger operator with a vr- 

periodic potential v G i7j^^(M). 

Then, for large enough n > N{v) the operator L has, in a disc of center and 

radius r„ = n/4, exactly two (counted with their algebraic multiplicity) periodic 

(for even n ), or antiperiodic (for odd n ) eigenvalues and X~ , and one Dirichlet 

eigenvalue fin- 
Let 

(6.5) A„ = |A+-A;^| + |A+-/i„|, n>N{v); 
then, for each sub-multiplicative weight uj and 

^=(Vt{n)), l^(n) = ^^, 

n 

we have 

(6.6) veH{n) {An)ef{n). 

Conversely, in the above notations, if uj = {uj{n))nez is a sub-multiplicative 
weight such that 

(p. Y^ logu;(n) 

(6.7) \ as n —>■ oo, 

n 

then 

(6.8) (A„) e e^{n) ^ ve H{n). 

If UJ is a sub -multiplicative weight of exponential type, i.e., 

(6.9) lim ^ ^ ^ > 

n— >oo fi 

then 

(6.10) (A„) G e^{n) ^ 3e > : v£ F(e^l"l). 

2. Throughout the paper and in Theorems 1221 and 1231 we consider three types 
of boundary conditions: Per^ and Dir in the form (a*) , (b*) and (c* = c) 
adjusted to the differential operators (jl.ip with singular potentials v G H~^. It 
is worth to observe that if v happens to be a regular potential, i.e., v G L'^{[0, vr]) 
(or even v G H°', a > —1/2) the boundary conditions (a*) and (5*) automatically 
become equivalent to the boundary conditions (a) and (6) as we used to write 
them in the regular case. Indeed (see the paragraph after (I2.67P ). we have 

(a* ) Per+ : y (vr) = y (0) , (y' - Qy) (vr) = {y' - Qy) (0) . 
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Therefore, with t> E L^, both the L^-function Q and the quasi-derivative u = 
y' ~ Qy continuous functions, so the two terms y' and Qy can be considered 
separately. Then the second condition in (a*) can be rewritten as 

(6.11) y'in) - y'iO) = Q(7r)y(7r) - g(0)2/(0). 
But, since Q is vr-periodic (see Proposition [1]) , 

(6.12) Q{7t) = Q(0), 

and with the first condition in (a*) the right side of ()6.1ip is Q(0)(y(7r)— y(0)) = 0. 
Therefore, (a*) comes to the form 
(a) y(7r) = y(0), y'{n) = y'(0). 

Of course, in the same way the condition (6*) automatically becomes equivalent 
to (6) ifvGH'', a> -1/2. 

A.Savchuk and A. Shkalikov checked ([27]), Theorem 1.5) which boundary 
conditions in terms of a function y and its quasi-derivative u = y' — Qy are 
regular by Birkhoff-Tamarkin. Not all of them are reduced to some canonical 
boundary conditions in the case of L^-potentials; the result could depend on the 
value of (5(0). For example, Dirichlet-Neumann be 

y(0)=0, (y'-Qy)(7r) =0 

would became 

2/(0) = 0, y'{TT) = Q(7r) • yin). 
Of course, one can adjust Q in advance by choosing (as it is done in [28] ) 



Q{x) = - v{t)dt if ?; G L^. 

J X 



But this choice is not good if Dirichlet-Neumann be is written with changed roles 
of the end points, i.e., 

(2/'-Qy)(0) =0, y{Ti) = Q. 

We want to restrict ourselves to such boundary conditions with v G that if 
by chance v G L'^ then the reduced boundary conditions do not depend on (5(0). 
We consider as good self-adjoint be only the following ones: 

Dir : y(0) = 0, y(7r) = 

and 

y(vr) = e'%{0) 

{y' - Qy){^) = e'\y' - Qym + Be'%{G), 
where 9 G [0, 27r) and B is real. 

Observations of this subsection are quite elementary but they would be im- 
portant if we would try to extend statements like Theorem [23] by finding other 
troikas of boundary conditions (and corresponding troikas of eigenvalues like 
{A+,A~,/i}) and using these spectral triangles and the decay rates of their di- 
ameters to characterize a smoothness of potentials v with a priori assumption 
V e (or even v G L^([0,7r]). 
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